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Abstract. We investigate the low-temperature phase diagram of the exactly solved 
su(4) two-leg spin ladder as a function of the rung coupling J± and magnetic field H 
by means of the thermodynamic Bethe Ansatz (TBA). In the absence of a magnetic 
field the model exhibits three quantum phases, while in the presence of a strong 
magnetic field there is no singlet ground state for ferromagnetic rung coupling. For 
antifcrromagnetic rung coupling, there is a gapped phase in the regime H < H c i, a 
fully polarized gapped phase for H > H c2 and a Luttinger liquid magnetic phase in the 
regime H c i < H < i/ c2 . The critical behaviour derived using the TBA is consistent 
with the existing experimental, numerical and perturbative results for the strong 
coupling ladder compounds. This includes the spin excitation gap and the critical fields 
H c i and -ff c2 , which are in excellent agreement with the experimental values for the 
known strong coupling ladder compounds (5IAP) 2 CuBr4-2H 2 0, Cu 2 (C5Hi 2 N 2 )2Cl4 
and (C5Hi 2 N) 2 CuBr 4 . In addition we predict the spin gap A J ± — | Jy for the weak 
coupling compounds with J± ~ J«, such as (VO) 2 P 2 Oy, and also show that the gap 
opens for arbitrary Jl/Jm. 



PACS numbers: 75.10.Jm,64.40.Cn 



TBA Phase diagram of the strong coupling ladder compounds 



2 



Recently there has been considerable theoretical and experimental interest in spin 
ladder systems. With the rapid progress presently being made in nano-engineering, 
many compounds with a ladder structure have been experimentally realized, such as 
SrCu 2 (B0 3 ) 2 , Cu 2 (C 5 H 12 N 2 ) 2 Cl 4 , (C 5 H 12 N) 2 CuBr 4 and KCuCl 3 [1. The existence 
of a spin gap, magnetization plateaus, superconductivity under hole doping, etc, are 
examples of some interesting physical properties that may be observed in experiments 
involving ladder compounds (see, e.g., Refs PP-[Z] and references therein). From the 
theoretical point of view, most of the results for ladder systems were initially obtained 
from studies of the standard Heisenberg ladder, which in contrast to its one-dimensional 
counterpart, cannot be solved exactly. Subsequently, other generalised ladder models 
have been proposed [H] and analysed through various numerical, approximate and exact 
approaches EDI CB EG] . 

On the other hand, although some exactly solved or integrable ladder models have 
been introduced jTSl UH UH| , none have been used to predict physical properties which 
could be compared directly with experimental data, such as the critical magnetic fields. 
In this context, the integrable spin ladder model based on the stt(4) algebra [13J appears 
to be a good candidate for this purpose, since its Hamiltonian consists of the standard 
Heisenberg ladder model with an extra biquadratic spin interaction term along the 
legs, the physical importance of which has been noted [S]. In the strong coupling 
limit, the contribution to the low-temperature physics from the biquadratic term is 
minimal, and as a consequence, the model exhibits similar critical behavior to the 
standard Heisenberg ladder. Therefore it is reasonable to expect that the integrable 
six (4) ladder model can well describe the low-temperature critical behavior of the strong 
coupling ladder materials. In addition, by properly minimizing the intrachain coupling 
in the integrable ladder Hamiltonian, the model may also be used to describe the weak 
coupling compounds. 

Here we show that this is in fact true in the strong coupling regime by investigating 
the quantum phase diagram of the integrable six (4) ladder, which can be tested by 
experiments. Our analytic expression for the gap, A = J± — 4Jj|/7, and the critical fields, 
fi B gH c i = A and /x#gif c2 = J± +4Ju/j, where 7 is a rescaling constant, can be applied 
in general to strong coupling ladder compounds with Heisenberg interactions, such as 
(5IAP) 2 CuBr 4 • 2H 2 0, Cu 2 (C 5 H 12 N 2 ) 2 Cl 4 , (C 5 H 12 N) 2 CuBr 4 and KCuCl 3 , by choosing 
7^4. For weak (J± ~ Jy) coupling compounds, such as (VO) 2 P 2 07, the choice 
7 ~ 8 determines a good fit for the gap PU E2] . In addition, in the presence of a strong 
magnetic field, we shown that the quantum phase diagram and the critical behavior 
predicted from the thermodynamic Bethe Ansatz (TBA) are in a good agreement with 
the experimental results for the above-mentioned compounds. We also show that the 
gap opens for an arbitrary value of J±/J\\, in accordance with the experimental results. 

I. The Model. We consider the phase diagram of the simplest integrable spin ladder 
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H = -H\ e( r 
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Here Sj and Tj are the standard spin-^ operators acting on site j of the upper and 
lower legs, respectively, Ju and J± are the intra (leg) and interchain (rung) couplings 
and h is the magnetic field. Throughout, L is the number of rungs and periodic 
boundary conditions are imposed. Essentially, the competition between the rung 
and leg couplings and the magnetic field h determines the physical properties and 
the critical behavior of the system. In order to facilitate the comparsion with real 
compounds, the intrachain part of this model (J2J) can be minimized through a rescaling 
constant 7. In comparison with the standard spin-| Heisenberg ladder [TJ [21 El El, 
the above Hamiltonian contains a four-spin interaction term, which minimizes the 
Haldane phase [Sj and causes a shift of the critical value of the rung coupling J± 
at which the model becomes massive. It is well established that this Hamiltonian is 
integrable and its leg part Hi eg is (up to a constant) simply the permutation operator 
corresponding to the su(4) algebra [13]. In addition, after the convenient change of basis, 
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TT>, |3> = 75 (I TO + I IT)) , |4) = I II), where the first 
state denotes the rung singlet and the three others the components of the triplet, the 
leg part remains of the same form while the rung term becomes diagonal. This rung 
term reduces the su(A) symmetry of H\ cg to su(3) © symmetry. Switching on the 
magnetic field breaks this symmetry further due to Zeeman splitting. This Hamiltonian 
can be diagonalized using the nested algebraic Bethe Ansatz (BA) with three levels. It is 
worth noting that for the ladder Hamiltonian ([TJ, the singlet rung state is energetically 
favoured for J± > 0, whereas the triplet rung state is favoured for J± < 0. On applying 
the magnetic field, component \2) of the triplet is energetically favoured. We will use 
these properties to our advantage by doing calculations with the choice of ordering for 
which the BA reference state is the closest to the true groundstate of the system. 

The underlying BA equations for Hamiltonian (JTJ) are well known JTSJ and consist 
of a set of three coupled equations depending on the flavors, v, u and w. Adopting the 
string conjecture ^3120] and taking the thermodynamic limit, the densities of the three 
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where pn (v), pn^ h (u) and pn^ h (w) denote the hole densities and 
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The symbol * denotes convolution and rjn (A) = 
1, 2, 3. The dressed energy el? plays the role 



of an excitation energy measured from the Fermi level. The driving matrix G depends 
on the choice of the reference state. Explicitly, for J± < 0, G = column(— ^-2ira n + 
nh, nh, —n(J± + h)) giving the free energy 

F(T, h) 



poo 00 (1) 

-h-T J2 a «( A ) ln ( X + e~^)dA. 
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On the other hand, for J± > 0, G = colum( — -2na n + n(J± — h),nh,nh), which leads 
to the form of the free energy (j7J) without the field term h. The TBA equations (j3*j) 
provide a clear physical picture of the groundstate, the elementary excitations, as well 
as the thermodynamic quantities, such as the free energy, magnetization, susceptibility, 
etc. Our results extend the earlier calculations on this model |13| I21j. 

II. Ferromagnetic rung coupling. In the low-temperature regime T — > 0, only the 
negative part of the dressed energies e"\ denoted by , contribute to the ground-state 
energy. The TBA equations (J3J) then become 
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g 3 _ a2 * e (3)- + ai * e (2)- ; ( 8 ) 

where g a , a = 1,2, 3, are the driving terms with respect to the basis order. In the 
regime J± < 0, the component | ft) of the triplet state is chosen as the reference 
state. The driving terms are given by g% = ——2iiai + h, g 2 = h and g% = —h — J±, 
respectively. Thus, in the absence of a magnetic field, the triplet is completely 
degenerate. The Fermi surface of the singlet is lifted as J± becomes more negative. If 
J_i_ is negative enough, the singlet rung state is not involved in the groundstate, namely 
e( 3 )(0) > 0, whereas the two other triplet Fermi seas still have their Fermi boundaries 
at infinity. In such a configuration, we may determine the critical point defining the 
transition from the su(4) phase into the su(3) phase by solving the TBA equations 
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(JSJ), with result J~ = — — In 3). At this critical point the free energy is given 

by F< ~°'°^ ~ — ^-(^(l) — V'Cl))) indicating a standard sw(3) phase. Here ^(n) is the 
digamma function. It is worth noting that the critical point J~ does not stabilize if an 
external magnetic field is applied. If the magnetic field is large enough, the ferromagnetic 
state |t|) becomes the true physical groundstate, i.e., there is a fully polarized gapped 
phase. It is found that for h > Hj = — , the state is fully-polarized provided that 

J_i_ < -. Therefore, in the ferromagnetic regime, the groundstate is in the critical 

su(3) phase. If the magnetic field is greater than H^, the groundstate is ferromagnetic 
with a magnetization plateau S z = 1. 

III. Strong antiferromagnetic regime. In the antiferromagnetic regime, J± > 0, 
the rung singlet state is the reference state. Thus the driving terms are given by 
gi = — ^27rai + J± — h and g 2 = g% = h, respectively. From the TBA equations (JHJ), if 
h = we immediately conclude that the triplet excitation is massive with the gap given 
by A = J_l — 4J||/7 for the regime J± > = — . Here is the critical point at which 
the quantum phase transition from the three branches of the Luttinger liquid phase to 
the dimerized u(l) phase occurs. To obtain good agreement with the experimental gap, 
we fix the rescaling constant 7 with the coupling constants remaining arbitrary. For 
the strong coupling compounds, e.g., (5IAP) 2 CuBr 4 • 2H 2 [5 , Cu 2 (C5H 12 N 2 ) 2 Cl 4 J3], 
(C 5 H 12 N) 2 CuBr4 J3], the experimental gap is well established as A « J_l — Ju, and as 
a consequence, we fix 7 « 4. On the other hand, for weak coupling compunds, e.g., 
(VO) 2 P 2 C>7 the choice 7 ~ 8 determines a good fit with the g U±. We 
stress that the purpose of introducing the rescaling constant is to minimize the effects 
of the biquadratic term, so that the model lies in the same Haldane phase as the pure 
Heisenberg ladder. 

IV. Magnetization plateau. The phase diagram of the antiferromagnetic spin ladders 
in the presence of a magnetic field is particularly interesting, because the critical 
points can be measured through critical magnetic fields. The appearance of quantized 
magnetization plateaus in the presence of a strong magnetic field is expected on general 
grounds PQ. From the TBA equations (JSJ) for antiferromagnetic rung coupling we observe 
that the magnetic field lifts the Fermi seas of and e^ 3 -*. If J± > J+, we can show 
that the two components of the triplet states, |3) and |4), do not become involved in the 
groundstate for a strong magnetic field. Basically, the magnetic field lifts the component 
1 2) of the triplet closer to the singlet groundstate such that they form a new effective 
spin-i state. Therefore, in a strong magnetic field the groundstate may be considered as 
a condensate of su(2) hard-core bosons. The gap can be deduced via the magnetic field h: 
the first critical field occurs at H c \ where g^sHd = A, i.e. the magnetic field closes the 
gap. The quantum phase transition from a gapped to a gappless Luttinger phase occurs. 
However, by continuing to increase the magnetic field h above the first critical field H&, 
the component |2) of the triplet becomes involved in the groundstate with a finite 
susceptibility. If the magnetic field is greater than the rung coupling, i.e. h > J±, the 
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state 1 2) becomes the lowest level. Therefore, it is reasonable to choose the basis order as 
(|2), |1), |3), |4)) T . Subsequently the driving terms are given by g 1 - 1 ' = —2irJ^ai — J± + h, 
g^ = J± and g( 3 ' = h. From the TBA, we see that the groundstate is a fully- 

4Jii 

polarized ferromagnetic state when the magnetic field is greater than H c2 = J± + —. 
Indeed, the critical field H c2 is in excellent agreement with the experimental data for 
the very strong coupling compound (5IAP) 2 CuBr 4 • 2H 2 (abbreviated as B5i2aT), 
and in a good agreement with the strong coupling compounds Cu 2 (C 5 H 12 N 2 ) 2 Cl4 
(abbreviated Cu(Hp)Cl) and (C 5 H 12 N) 2 CuBr 4 (abbreviated BPCN) @] (see Table 
1). On the other hand, the precise structure of the compound KCuCl 3 is not clear pQ. 
It is believed to exhibit a double-chain structure [0] with a gap A « 35K identified 
via the best fitting in the susceptibility curve through the Troyer formula [22]- The 
coupling constants are determined as J± = 4Jy, J\\ = 12. 3K, Jdi ag = jHj. However, 
high-field measurements indicate the gap A rs 31. IK [Zj. Our TBA result gives poor 
agreement with the experimental result for this type of ladder compound (see Table 1). 
This suggests that the compound may exhibit a double-chain structure with additional 
diagonal interaction. For these double-chain structure ladders, such as KCuCl 3 , TlCuCl 3 
etc., the leg couplings appear to be very large, resulting in discrepancy with the critical 
fields derived from the TBA method. 

After a similar calculation, we obtain the magnetization S z ~ 4Qi(l — 
2Qi/tt)/ti in the vicinity of the critical field H c i, with the Fermi boundary Q\ « 
v/(/i — H c i) j (H c i — 5h) . For a very strong magnetic field such that H c2 — h « 1 
the free energy is 



*(f--*)y (9) 

and the susceptibility k ta -^=j= (H C 2 — h)~^, which indicates the nature of 
the singular behavior in a phase transition from a gapless to a ferromagnetic 
phase. The magnetization is given by S z ~ 1 — 4Q 2 (1 — 2<5 2 /7r)/7r, where Q 2 ~ 
a/ (H c2 — h)/ (5h — H C 2). The fact that the magnetization depends on the square root of 
the field in the vicinity of the critical fields is consistent with other theoretical fTJ ^] 
and numerical results jHj- The magnetization increases almost linearly between the 
critical field H c i and H c2 . The groundstate is ferromagnetic above H c2 with the gap 
A = fig(H-H c2 ). 

Numerical solution of the TBA equations gives a reasonable magnetization curve 
(see Fig. 1) which passes through an inflection point midway between if cl and H c2 . This 
inflection point is clearly visible in experimental curves, e.g. for (5IAP) 2 CuBr4-2H 2 5J, 
Cu 2 (C 5 H 12 N 2 ) 2 Ci4 [3] and (C 5 H 12 N) 2 CuBr 4 jlj. The physical meaning of the inflection 
point is that the probabilities of the singlet and the triplet states |2) in the groundstate 
are equal. It suggests an ordered dimer state close to half-filling [2H]- Therefore, in the 
strong-coupling regime, the one point correlation function (Sj-Tj) = — | lies in a gapped 
singlet groundstate, which indicates an ordered dimer phase, while (Sj ■ Tj) — | is in 
the fully-polarized ferromagnetic phase. However, in a Luttinger liquid phase, we find 
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Compounds 


9 


J± 


J|, 


7 




F c2 (exp) 


# c i(TBA) 


# c2 (TBA) 


B5i2aT 


2.1 


13K 


1.15K 


4 


8.4T 


WAT 


8.3T 


10.03T 


Cu(Hp)Cl 


2.03 


13.2K 


2.5K 


4 


7.5T 


13. 2T 


7.84T 


11.51T 


BPCB 


2.13 


13.3K 


3.8K 


4 


6.6T 


14.6T 


6.6T 


11.95T 


KCuCl 3 


2.05 


49.2K 


12.3K 


2.68 


22.4T 


« 60T 


22.4T 


49T 



Table 1. Comparison between the experimental values for the critical points H c \ and 
H C 2 for strong coupling ladder compounds and the TBA results obtained from the 
su(A) integrable model. 



(Sj : ■ Tj) = — | + S z . The magnetic field increases the one point correlation function. 

We also notice that our results for the gap, A = J± — AJ^/j, and the critical 
field, H C 2 = J± + 4J11/7, coincide for 7 = 4 with the first-order perturbation theory 
results obtained for strong coupling [IT] . However, their higher-order terms lead to poor 
agreement with the experimental results. It is apparent that the rescaling constant 7 
causes a shift in the critical point. This can be seen from the values of 

H c2 /A = 1 + 2/(^-1), (10) 

which are plotted in Fig. 2. The larger the ratio of y-, the closer the two critical 
points are. This means that the critical points H c \ and if c2 cannot be distinguished for 
a very large energy gap. Once the gap is closed by an external field, the groundstate 
immediately becomes fully-polarized. This is evident in the strong coupling compound 
(5IAP) 2 CuBr 4 • 2H 2 5|. Here the gap opens only if J±/J\\ > 4/7, with 7 arbitrary. 
And therefore the gap opens for arbitrary J±/J\\. 



Finally, we show the phase diagram in the presence of a magnetic field in Fig. 3. In 
the ferromagnetic rung coupling regime, the fully-polarized ferromagnetic state lies in 
the region h > \J±\ and h > 4J11/7, whereas the su(3) Luttinger magnetic phase is in the 
region h < \J±\ and left of the boundary between the su(3) and sw(4) phases. The sit(4) 
phase is in the region h < 4Jy/7 and right of this boundary. In the antiferromagnetic 
rung coupling regime, the singlet rung state lies in h < J± — 4J||/7 whereas the 
ferromagnetic fully-polarized state is in the region h > J± +4J11/7. The su(2) magnetic 
phase remains in the region h > J± — 4J11/7, h < J± + AJn/^f and J± > 4J11/7. The 
sw(4) magnetic phase lies in the region h < J± + 4Jy/7 and < J± < 4:J\\/j. 

To conclude, we have studied the phase diagram of the integrable sw(4) spin ladder 
model (1) by means of the TBA. In particular, the critical behavior at the critical points 
H c i and H c2 was derived. In the presence of a strong magnetic field, the phase diagram 
is in a good agreement with the experimental observations for the strong coupling 
compounds (5IAP) 2 CuBr 4 • 2H 2 5J, Cu 2 (C 5 H 12 N 2 ) 2 Cl 4 and (C 5 H 12 N) 2 CuBr 4 i). 
We have also predicted the spin gap A w J ± — ~J|| for the weak coupling compounds 
with J± ~ J11, such as (V '0) 2 P 2 07 and also shown that the gap opens for arbitrary 
J±/J\\- 
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Figure 1. Magnetization S z versus magnetic field H — Hugh obtained from the TBA 
equations for the values J± = 13K, Ju = 1.15K and 7 = 4 for the strong coupling 
compound (5IAP)2CuBr4 • 2H2O 5 . At the inflection point h = J± the magnetization 
is 0.5. The curve is in excellent agreement with the experimental result [S]- 






1 3 5 7 9 11 13 15 



Figure 2. The critical point H c2 /A as a function of the ratio J±/J\\ for different 
values of the rescaling parameter 7. Also shown is the perturbation theory result. 
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Figure 3. The magnetic phase diagram of the two-leg su(A) ladder. In the 
antiferromagnetic regime the thick lines are h = J± — 4J11/7, h — Jj_ + 4J||/7 and 
the dash line is h — J±. In the ferromagnetic regime the thick lines are h = — J± 
and h = 4J||/7. The dash-dot line is an approximate boundary between the su(4) and 
su(3) phases. 
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